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graphite with Bernal stacking structure is a semiconductor with a narrow band gap 
Eg - 40 meV. 
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1. Introduction 

The study of the transport properties of graphite has been object of discussion for 
more that 60 years [H [21 [3l HI O |6] . The work of Slonczewski and Weiss |lj as well as 
that of McClure [5l |6] - both a continuation of the work of Wallace [1] on the band 
structure of graphene - introduced several (free) coupling parameters within a tight 
binding calculation and the k.p. method to obtain the graphite band structure. The 
parameters were then fixed using as reference transport results obtained in bulk graphite 
samples with carrier density per graphene layer ?7,(T — ?■ 0) = rio ^ 10^*^ cm~^. The 
apparent existence of a relatively large carrier density seemed to be compatible with the 
relatively small and metalliclike resistivity measured in several bulk graphite samples 
[7]. All these early results obtained for the graphite structure were resumed later within 
the so-called two-bands model for graphite developed by Kelly [7]. After that, many 
band structures calculations using the local density as well as others approximations 
[8] reached the same conclusions and claimed to fit well the experimental results, see 
e.g. Refs. [9],[T0]. Note, however, that all these calculations depend on parameters that 
are being fixed to fit experimental data. If those experimental results do not reflect the 
intrinsic properties of ideal graphite, clearly those calculations as well as the obtained 
parameters cannot be taken as intrinsic of the graphite structure. 

A different approach to obtain the binding energy between the planes of graphite has 
been proposed using Lennard- Jones potential and assuming that the binding between 
graphene layers is due to van der Waals forces [HI [12], [13]. These forces represent a 
very weak interaction and therefore it is rather clear that this approach will not give 
a similar band structure for graphite as the one proposed with large coupling between 
the graphene planes. From the experimental point of view we may now doubt that the 
large amount of the reported data, e.g. n{T), do reflect ideal graphite. The exhaustive 
experience accumulated in gapless or narrow gap semiconductors [11] indicates already 
that care should be taken with the measured carrier density because it can easily be 
influenced by impurities and/or defects in the graphite/graphene lattice [T5l [T6]. 

We believe that the multigrain morphology of real oriented graphite samples has 
not been taken into account appropriately in the literature, especially its influence to the 
carrier density and other transport phenomena. For example, electron back scattering 
diffraction (EBSD) reveals typical size of the grains of a few microns in the a, b plane 
within HOPG samples, see Ref . [17] , a size that limits the carrier mean free path at low 
temperatures [18]. Moreover, the contribution to the transport properties of internal 
two-dimensional interfaces found recently between single crystalline regions inside the 
HOPG samples [I9] was apparently never taken into account in the existence literature 
of graphite. Why these interfaces can be of extreme importance in graphite? As the 
results from various semiconductors show (e.g. n-Ge bicrystals [201 [21], p-InSb [22] as 
well as in Hgi^^jCdj-Te grain boundaries [23]) internal interfaces lead to the formation 
of confined quasi-two dimensional carrier systems with uq ~ 10^^ . . . 10^'^ cm~^ and 
clear signs for the Quantum Hall effect [21]. Recently done studies demonstrated the 
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large sensitivity of the resistivity of graphite samples to the internal interfaces that 
exist between crystalline regions of ~ 30 . . . 100 nm thickness, a few microns long and 
running mostly parallel to the graphene planes [T9]. These results indicate that the 
earlier reported values of n{T) as well as the metalliclike behavior of the resistivity p{T) 
are not intrinsic of the graphite Bernal structure but are due to a large extent to the 
contribution of internal interfaces and defects. 

The aim of this work is to propose a simple model to explain the experimental 
longitudinal resistance data obtained in different oriented graphite samples of different 
thickness and areas. The rather complicated behavior of the longitudinal resistivity 
can be explained assuming the parallel contribution of regions with semiconducting 
graphene layers and the ones from the interfaces between them, as transmission electron 
microscopy studies revealed [T9]. 

2. Sample characteristics 

All the graphite samples were taken from the same Bernal-type highly oriented pyrolytic 
graphite (HOPG) sample from Advanced Ceramic of high purity (see, e.g., Ref. |24j ) 
and with a rocking curve width of (0.35 ±0.1)°. As shown by Bernal [25]. this structure 
has the usual ABABABA. . . stacking of the single graphene layers inside graphite. This 
is the stable graphite structure, which is obtained by suitable annealing treatments 
eliminating also possible rhombohedral modification that results from deformation of 
the original hexagonal structure. The internal structure of the used samples is shown 
in a transmission electron microscope picture in Fig. [H As shown in Ref. [I9] this 
picture reveals single crystalline regions of graphene layers of thickness between 30 nm 
and ~ 100 nm. These regions are clearly recognized in Fig. [T] through the different 
gray scales and have sightly different orientation between each other, e.g. different angle 
misalignments. We note that a rotation up to 30 degrees between the graphene layers 
from neighboring graphite regions has been seen by HRTEM in few layers graphene 
sheets in Ref. [26] . 

Important to understand the measured behavior of the resistance of graphite 
samples is the existence of well defined interfaces between the single crystalline 
regions. As we mentioned in the introduction, interfaces between crystalline regions in 
semiconductors with different orientations for example [20 1 [2T 1 [22 l [23]. lead to confined 
quasi-two dimensional carrier systems with much larger carrier density than the bulk 
matrix. We assume therefore that these interfaces running parallel to the graphene 
layers of the graphite structure are the origin of the metallic like resistivity as well as 
for the apparent large carrier density measured in bulk samples, i.e. hq > 10^*^ cm~^, an 
assumption that is supported by the change of the absolute resistivity with thickness 

Taking into account the internal microstructure of the graphite samples it is clear 
that it is necessary to measure the resistance of samples of small enough thickness in 
order to get the intrinsic transport of the graphite structure with its weak coupled 
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Figure 1. Transmission electron microscope picture of the internal microstructure of a 
highly oriented pyrolytic graphite sample as shown in Ref. |19j and similar to those used 
in this study. The dashed red lines indicate some of the interface regions between the 
single crystalline graphite parts. These interfaces run parallel to the graphene layers 
up to several micrometers till a larger defect is encountered (see the defective region in 
the middle of the picture) . Regions with different gray colors indicate slightly different 
orientation of the graphite structure. From measurements with EBSD we know that 
the usual size of these single crystalline regions in the a, 6 plane is less than 10 /.jm 

mm. 



graphene layers. The micrometer small samples and with thickness below 100 nm 
were prepared by a rubbing method [T9]. The resistivity data were obtained with four 
contacts and checking the ohmic behavior at room temperature. The deposited Pd/Au 
contacts on the micrometer large graphite flakes were prepared by electron lithography 
and micro- Raman characterization was used to check for the sample quality [19], as 
the inset in Fig. [2] demonstrates. Following the characterization work of Ref. |27], the 
Raman spectra confirm the similar stacking between the bulk and thin graphite samples. 
Micrometer large mean free path at 300 K and low carrier density [28| [29] are a further 
proof of the high quality of the thin graphite flakes discussed in this work. 

3. Temperature dependence of the carrier density 

We start with the temperature dependence of the carrier density n{T) obtained from 
a 40 nm thin and several micrometers large graphite flake sample. It is important 
to recognize that due to the large mean free path of the carriers I in the graphene 
layers within the graphite structure, ballistic not diffusive transport has to be taken 
into account [T8l [29] . To obtain the mean free path or the intrinsic carrier density of 
the graphene layers inside graphite one cannot use straightforwardly the Boltzmann- 
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Drude approach to interpret the longitudinal and transverse resistances. Due to the 
large mean free path of the carriers in the weakly coupled graphene layers within 
the graphite structure it is possible to use experimental methods where the ballistic 
transport is clearly revealed. One possible experimental method is the constriction 
method based on the measurement of the longitudinal resistance i? as a function of 
the width W oi a, constriction located between the voltage electrodes. When i ^ W 
the ballistic contribution overwhelms the diffusive ones and allows us to obtain i{T) 
and n{T) without the need of free parameters or arbitrary assumptions [IHl [28]. Other 
experimental method is the measurement of the length dependence of the resistance. If 
ballistic transport is important then a finite resistance is extrapolated at zero sample 
length, a value that can be used to obtain directly the mean free path without free 
parameters or arbitrary assumptions [29]. Both independent methods provide similar 
large carrier mean free path as well as much lower carrier density as those found in the 
literature for bulk graphite samples. 

Figure Elshows n{T) for a thin graphite sample; it follows an exponential dependence 
with an activation energy or energy gap Eg ~ 46 meV. Due to the small gap and from 
Fig. [2] it is not obvious that the experimental curve follows an exponential behavior and 
not a linear one at temperatures above the saturation region. A simple way to distinguish 
that, it is calculating the temperature derivative. In case of a linear temperature 
behavior we expect a saturating upper derivative value at high temperatures but a 
shallow maximum for the exponential function at intermediate temperatures. The 
bottom right inset shows the calculated derivative for both curves, the experimental 
and the exponential function fit. It is clear that the exponential function provides the 
correct temperature dependence of the carrier density. The shallow minimum in R{T) 
(main panel in Fig. [2]) at T < 50 K is an artifact due to interfaces, at least that one 
between the substrate and sample and/or the sample free surface. The finite value of 
n{T — > 0) can be due to this interface contribution or due to lattice defectsfTH] including 
a very small amount of impurities like hydrogen; note that n = 10^ cm~^ would mean 
of the order of a single hydrogen atom or C-vacancy in 1 /im^ graphene area. 

One may be surprised to get such low carrier density in comparison with usual values 
reported in the literature for graphite, e.g. Uq > 10^° cm~^ [3 El [10]. We note first that 
the measured graphite sample has a thickness less than 50 nm and therefore much less 
contribution of interfaces. Second, the measured carrier density of graphite reported 
in literature was obtained mostly for bulk samples with an unknown concentration of 
interfaces as well as defects. Further independent evidence that supports an intrinsic 
low carrier density in graphite is given by the vanishing of the Shubnikov-de Haas (SdH) 
oscillations in the magnetoresistance the thinner the graphite sample, an experimental 
fact already reported in 2001 that to our knowledge remained without explanation 
[30] . Note that older publications reported SdH oscillations of very large amplitude, 
e.g. Ref. [31], in comparison with the rather weak amplitudes, if at all observed, 
found nowadays in better quality samples or thin enough graphite samples. A clear 
demonstration that these SdH oscillations are actually non-intrinsic of ideal graphite 
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Figure 2. Carrier density per graphene layer obtained using the constriction method 
[l8] for a graphite sample of size 9 x 3 x 0.040 ^td? . The continuous line is a fit to 
the data and follows 7i(r)[10^cm"2] ^ 2.7 + 12.4 exp(-540/2T[K]). We estimate a 
~ 30% error in the absolute value of the carrier density, mainly due sample geometry 
errors as well as in the constriction widths. The upper left inset shows the Raman 
(514 nm) spectra of: (o) bulk graphite (Fig. 2(a)) and a (red line) multigraphene 
sample (similar to that of Fig. 3(a)). The absence of a D-peak indicates the absence 
of a significant number of defects. The bottom right inset shows the temperature 
dependent derivative of the carrier density (circles) and of the fitting curve shown in 
the main panel. Note the maximum in the derivative at 125 K present in both 
derivatives. The horizontal straight line is only a guide to the eye. 



is given by their appearance in thin graphite samples after introducing defects by 
irradiation as shown in Refs. [T^ [32]. 

4. Temperature dependence of the resistance of different graphite samples 

Following the results from Ref. [19] as well as the semiconducting behavior of n{T), see 
Figl2l obtained for a thin graphite sample, we assume therefore that the graphene layers 
inside each graphite sample are semiconducting and altogether and between the voltage 
electrodes their signal depends on an effective resistance of the type: 



The prefactor a(T) depends basically on the mobility, i.e. the mean free path and on 
details of the carriers band structure (e.g. effective mass). If we take into account in 
a(T) the temperature dependence of the mean free path recently obtained for similar 



samples, i.e. i{T) ~ ((3)-^ + (6.4 x lO^T^)-^)-^ {i in /zm and T in K), see Ref. [2l 




the absolute value of the energy gap Eg obtained from the fits to the data remains 




(1) 
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the same within the confidence hmits of 35%. For simphcity we will take a(T) as a 
temperature independent parameter as well as the energy gap Eg. The absolute value 
of the prefactor a (as well as of the other prefactors) will change from sample to sample; 
note that we estimate the resistance of a given sample not its resistivity. 

For samples with thickness larger than ~ 50 nm and of several micrometers length 
there is a larger probability to have interfaces, which signals will be picked up by the 
voltage contact electrodes (usually several micrometers apart). Therefore, in parallel to 
RsiT) we simulate the contribution from the interfaces through the resistance: 

R,{T) = Ro + R,T + R2 eM-Ea/keT) , (2) 

where the coefficients R2 as well as the activation energy Ea are free parameters. The 
temperature independent term in Eq. (jj]) represents the residual resistance measurable 
at low enough temperatures. The unusual thermally activated term as well as the, in 
general weak, linear one will be discussed at the end. The total resistance Rt{T) is 
given by the parallel contribution of Rs and Ri as Rt{T) = [Rs(T)~^ + _Rj(T)~^]~^. 
Clearly, by changing the parameters one can obtain all types of behavior for R{T). We 
will see that a consistent description of the data can be indeed reached and that the 
main free parameter Eg is similar for all samples. 

Figures [3] and m show the normalized resistance vs. temperature and the fits with the 
parallel resistance model described above of six different samples with different weights 
between the intrinsic semiconducting and the non-intrinsic interface contribution. That 
is why the pre-factors of the main terms in Rt{T) must change from sample to sample, 
i.e. Ri/ Ro, R2/ Rq and o/Rq, because we do not fit the resistivity but the resistance. 
From all the obtained data and within the confidence limits of the fitting we obtain an 
energy gap Eg = AO ± 15 meV, independently of the sample geometry, i.e. a similar 
energy gap is obtained for very thin as well as thick graphite samples if the contribution 
of the interfaces does not short circuit completely that from the graphite crystalline 
regions. This fact indicates that we are not dealing here with a special graphite structure 
found only in a few tens of nanometers thick samples. This is the main result of this 
study. 

We stress that the interfaces as well as the single crystalline regions are restricted in 
thickness as well as in the a, b plane parallel to the graphene layers. This fact supports 
the use of the simple parallel resistance model. For thick enough graphite samples 
(thickness t > 1 fim) the contribution of the interfaces overwhelms and the measured 
resistance will be mainly given by Eq. (j2]), as the results for bulk HOPG sample shows, 
see Fig. ^a). The exponential term with an activation energy Ea/ks = (48 ± 2)K is 
of the same order as reported in Ref. |33] for similar samples. To fit the data shown in 
Fig. El^a) a relatively small linear-temperature term contributes at high temperatures. 

In general the semiconducting behavior starts to be clearly visible in thinner samples 
due to the decrease in the amount of interfaces; obviously there is no sharp thickness 
threshold because the properties of the interfaces as well as the defects and carrier 
concentrations vary from sample to sample, affecting the relative weight between the 
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Figure 3. Normalized resistance R/Ra vs. temperature for three HOPG samples 
with thickness, length (between voltage electrodes), width and Rq: (a) ~ 20 iim, 
2 mm, 1 mm, 0.003 fi; (b) ~ 10 /im, 1 mm, 1 mm, 0.013 fl; (c) ~ 50 nm, 
~ 3 yum, ^ 3 /im, 15 fl. The (red) lines through the experimental data are 
obtained from the the parallel contributions given by Eqs. ([T]) and ([2]): (a) Ri/Ro = 
1 + 2.2 X 10-3r[K] + 5.2exp(-48/T[K]); (b) i?.,/i?o = 1 + 2.2 cxp(-33/r[K]) and 
Rs{T)/Ro = 3.4exp(662/2r[K]); (c) Ri{T)/Ro = l+2xlG-3T[K]+0.7cxp(-38/T[K]) 
and Rs{T)/Ro = 2.35 exp(340 /2T [K] ) . Note that the curves obtained are very sensitive 
to small changes of parameters. For example in (c), the continuous curve just below 
follows R.,{T)/Ro = 1 + 2.41 x 10-3r[K] + 0.7exp(-38/r[K]) and Rs{T)/Ro = 
2.35exp(260/2r[K]). Whereas the dashed curve above follows Ri{T)/R^ = 1 + 2.17 x 
10-3r[K] +0.7exp(-38/T[K]) and Rs{T)/Ro = 2.35 exp(400/2T[K]). The inset in (c) 
blows out the high temperature part where one recognizes the slight increase of the 
resistance at T > 200 K, an increase that can be obtained changing slightly the used 
parameters as the lower curve shows. 




Temperature T{K) Temperature T(K) Temperature T(K) 
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Figure 4. Similar to Fig. [3] but for samples with geometry and Rq: (a) ~ 13 nm, 
14 ^m , 10 /im, 490 il; (b) ~ 20 nm, 5 /im ,10 /im, 32 rj; (c) ~ 37 nm, 27 /im, 
6 fim, 69 il. The lines through the experimental data follow: (a) Ri{T)/Ro = 
1 - 6.2 X 10-''r[K] and i?s(T)/i?o = 1.18 exp(480/2T[K]); (b) R,(T)/Ra = 1 + 
0.096 exp(-24/T[K]) and R,{T)/Ro = 0.82 exp(350/2r[K]). (c) R,{T)/Ro = 1 + 
0.0044T[K] + 0.47cxp(-12/r[K]) and Rs{T)/Ra = 0.73 exp(400/2r[K]. 



two parallel contributions. The results shown in Figs. Et^b) and (c) are examples where 
the two contributions Ri and Rs compete. For samples thinner than ~ 40 nm the 
semiconducting behavior starts to overwhelm that from the interfaces, see Fig. HI 

Note that this semiconducting behavior is observed for a range of about five orders 
of magnitude in sample area and that the width or length of the samples overwhelms 
the intrinsic Fermi wavelength of the graphene planes inside the sample [28] Ai? < 1 /im. 
Therefore, the semiconducting behavior is not due to a quantum confinement shift, as 
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observed in Bi nanowires[M] of small diameter d < Xp. We stress that our results 
should not be interpreted as a metallic to semiconducting transition as a function of 
thickness, but only the weight between the two parallel contributions changes upon 
sample thickness. Because at least the interfaces sample-substrate and sample-surface 
remain, as for the sample of Fig. Ufa), both parallel contributions are always present. 
We note that the T— dependence of the c— axis resistivity of graphite with pressure 
[35] shows similar characteristics as the ones for the resistivity in plane shown here. 
The c— axis resistivity behavior was basically interpreted with a tunnel barrier for the 
electrons between planes and crystallites. 

5. Discussion 

We would like to discuss the question of the existence of an energy gap taking into 
account results from literature and whether this intrinsic semiconducting behavior has 
been observed in the past, (a) As we discussed above and taking into account the 
contribution of the interfaces bulk graphite samples have, it is clear that the behavior 
obtained in transport measurements of those samples is not intrinsic. Therefore we 
cannot consider those experimental results as a proof for the apparent metallicity of 
graphite. 

(b) One would expect that infrared (IR) measurements show somehow the 
semiconducting behavior. However, those measurements, as for example in Ref. [36], 
were done in bulk samples (several mm^ area and several hundreds of micrometers thick) 
that show clear SdH oscillations [37] with the usual carrier density above 10^° cm~^ per 
graphene layer and therefore do not represent intrinsic graphite. Note that the internal 
interfaces, see Fig. [T], and due to their large carrier density in comparison with the 
crystalline regions, are those regions that provide a shielding to the electromagnetic 
wave used in IR measurements. In general we remark that spectroscopy methods are 
not completely reliable when one wants to resolve a gap of 40 meV in samples that 
are inhomogeneous for the reasons discussed above. Therefore, we believe that the 
intrinsic properties of the single crystalline graphite regions are not actually seen by 
this technique. Also other experimental techniques, see e.g. Ref. [38], would have 
problems to recognize an electronic band structure with an energy gap of the order of 
40 meV if the carrier density is much larger than the intrinsic one of ideal, defect free 
graphite. 

(c) We may ask about the results on graphite using angle-resolved photoemission 
spectroscopy (ARPES) published several times in the past. There are three points 
we should take into account about ARPES results. Firstly, this technique is surface 
sensitive, although we may expect to see to some extent a vestige of the bulk band edge. 
Second, different samples may provide different results. A clear correlation between 
sample quality and other of its characteristics with the ARPES results is difficult because 
one not not always find the necessary information in the reported studies. Finally the 
scanned sample area and the energy resolution of the usual spectroscopy systems. We 
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know from other techniques hke electron back scattering diffraction (EBSD) [17j or 
electron force microscopy (EFM) |39] that homogeneous regions of graphite appear to 
be in the region of 10 /xm (parallel to the planes) or smaller. Regarding the energy 
resolution we find, for example, in the ARPES work of Ref. [IQ] a resolution of 15 meV to 
65 meV, in Ref. it is from 15 meV to 40 meV and in Ref. |12] 40 meV. This resolution 
appears to be too low to clearly resolve an energy gap of the order of 30 to 50 meV 
in bulk graphite. A better energy resolution between 4 to 15 meV has been achieved 
in Refs. |l3l |9] using kish and natural graphite samples, respectively. Interestingly, in 
Ref. [l3] an energy gap of 25 meV between the vr and vr* bands at the K(H) point 
was reported whereas in Ref. [9] an energy gap of 37 meV has been inferred from the 
band fits near the H point, although the used sample had apparently an average carrier 
density of ~ 6 x lO^^cm"^ (at 25 K). Clearly, higher energy resolution, smaller scanned 
areas and clearer sample quality aspects are necessary in future experiments using this 
technique. 

(d) In general we note that techniques that are sensitive to the surface do not 
necessarily provide the intrinsic properties of graphite. Apart from the extra doping 
graphite surfaces may have due to defects and added atoms or molecules, we note that 
the surface itself may have a different electronic structure as the graphene layers inside 
the Bernal structure due to the total absence of coupling with a graphene layer above 
it. Therefore, STM experiments on the graphite surface do not necessarily provide the 
intrinsic behavior of graphite with certainty. However, there are STM experiments that 
suggest the existence of an energy gap of the same magnitude as the one obtained in this 
work. This energy gap appeared in a few graphene layers [H] from "unidentified states" 
or levels that at zero magnetic field do not coalesce, a result supported by Hubbard 
calculations with a constant potential |15]. Clearly, if a few graphene layers shows an 
energy gap, there seems to be no simple reason to rule out that graphite should not 
have it. Note also that in Ref. other electronic contributions, partially due to Dirac 
fermions, have been obtained at the surface. This fact suggests that the resistance of 
this interface (or between the sample and substrate) will provide a different conduction 
path in parallel to the semiconducting one from the internal graphene layers. This gives 
a simple explanation for the saturation in resistance measured at low temperatures. 

(e) In order to compare the results of the temperature dependence of the resistance 
R{T) obtained in single or a few layers graphene samples with our results one should 
compare the one obtained at the lowest carrier density possible. This has been done in, 
e.g., Ref. [l6] where one recognizes that at the lowest carrier density reached near the 
neutrality point (former Dirac point) the resistance does not behave metallic but rather 
semiconducting. Increasing the carrier density only slightly with the bias voltage, the 
resistance starts to show signs of metallicity[l6] . Further semiconducting behavior in 
R{T) has been seen in Refs. |17j and ^8] . 

(f ) And at last but not least we would like to point out the recently done optical 
pump-probe spectroscopy experiments on 30 nm thick graphite flake [19], a sample 
geometry relevant for comparison with our study. Even for excitation densities 10 times 
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larger than the one of our samples, these experiments suggest a renormahzation of the 
band gap by ~ 30 meV, which is nearly constant during the first picosecond. The authors 
in Ref. [19] conclude that carrier equilibration in graphite is similar to semiconductors 
with a nonzero band gap. 

The interaction between graphene planes in graphite is van der Waals type, which 
gives a binding energy of < 400 K. That means that each graphene layer should be 
very little affected and therefore the electron electron interaction may be also active in 
graphite as also in ideal single graphene as well. This cannot be excluded because the 
data obtained in the literature are for large n and therefore any gap of the magnitude 
we obtain here will not be so easily measurable. Electron interactions are large and 
for a small enough carrier density the expected screening will be very weak promoting 
therefore the existence of an energy gap. This is what it is observed in Monte Carlo 
simulations for the unscreened Coulomb interaction in graphene with different Dirac 
flavors [50]. 

Before concluding we would like to discuss the possible origin of the exponential 
function in Eq. ([2]) because this is not the usual one expects for metals or semimetals 
and cannot be understood within the usual electron-phonon interaction mechanisms also 
in two dimensions. We note that this function has been already used to describe the 
increasing resistance of bulk graphite samples with temperature and speculated to be 
related to some superconducting- like behavior in graphite [33]. We note further that a 
similar dependence has been observed in granular Al-Ge[5T], which shows for a particular 
Al concentration a semiconductor-superconductor transition practically identical to that 
in Fig.|4](c). The thermally activated behavior in Eq. ([2]) can be understood on the basis 
of the Langer-Ambegaokar-McCumber-Halperin (LAMH) model [521 [53] that applies to 
narrow superconducting channels in which thermal fluctuations can cause phase slips. 
The value of the activation energy Ea obtained from the fits depends on the measured 
sample and it is between 10 K and 40 K for a large amount of measured samples. 

Although we do not claim that this exponential function is due to granular 
superconductivity embedded within the interfaces shown in Fig. [H this speculation 
matches however a series of experimental hints obtained recently. We would like to 
mention that the possible existence of high-temperature superconductivity in non- 
percolative regions within the internal interfaces found in highly oriented pyrolytic 
graphite has been proposed recently [19] . The experimental hints are quite diverse 
and all of them point in the same direction: (1) Huge magnetic field driven metal- 
insulator transition in bulk graphite [5H [55] . This transition vanishes for thin 
enough samples because of the lack of the internal interfaces. (2) Superconducting- 
like hysteresis loops in the magnetization of pure HOPG samples [561 ISZj- (3) 
Anomalous hysteresis in the magnetoresistance of mesoscopic graphite samples, similar 
to those found in granular superconductors [Ml [59] . (4) Quantum oscillations in the 
magnetoresistance in multigraphene samples based on Andreev-scattering mechanism 
[58] . (5) Superconducting-like magnetic field driven transition in the resistance of an 
internal interface of a graphite flake [60]. (6) Finally, Josephson-like I /V characteristic 
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curves found in thin graphite lamellas [HI]. 

Since most of the theoretical predictions emphasize that superconductivity should 
be possible under the premise that the carrier density per graphene area n > 10^^ cm~^ 
in order to reach > 1 K, it is then appealing that regions at the internal interfaces 
within the graphite structure may have superconducting patches due to its much higher 
carrier density. We note that interfaces in some semimetals materials like Bi can be 
superconducting [62], |63] . 

The weak linear term in Eq. ([2]) provides a typical metallic dependence to the 
resistance, which may come from the metallic regions within the interfaces. This small 
term is not always necessary to get a reasonable fit to the data and can be even negative 
(or may follow a weak variable range hopping dependence) as appears to be the case 
for the sample-substrate and/or sample-surface interfaces, see Fig. ID^a). Even in the 
case of total absence of internal interfaces, as for example for thin enough samples, it 
is clear that the semiconducting behavior of the intrinsic graphite regions cannot be 
observed at low enough temperatures due to the parallel contributions of the surface 
and/or substrate/sample. 

Concluding, the temperature dependence of the electrical resistance of oriented 
graphite samples can be quantitatively understood assuming the parallel contribution 
of semiconducting and normal metallic (and/or granular superconducting) regions. The 
semiconducting contribution clearly indicates that Bernal graphite is a semiconductor 
with a narrow gap of the order of 40 meV. We speculate that the existence of narrow 
band gap we observe in graphite may be applicable to other semimetals. In other words 
we doubt on the existence of strictly zero band gap in nature. 
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